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STABILITY FOR SOME INVERSE PROBLEMS FOR TRANSPORT 

EQUATIONS 

FIKRET GOLGELEYENt AND MASAHIRO YAMAMOTO* 

Abstract. In this article, we consider inverse problems of determining a source term and a 
coefficient of a first-order partial differential equation and prove conditional stability estimates with 
minimum boundary observation data and relaxed condition on the principal part. 
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1. Introduction and main results. Let fl C M" be a bounded domain with 
smooth boundary dfl and let I'ix) be the unit outward normal vector to dil. Let us 
consider 

dty{x,t) + H{x) ■ Vy{x,t) + V{x)y{x,t) = f{x)R{x,t), xGV,,0<t<T (1-1) 
and 


j/(a;,0) = 0, X £ n. (1-2) 

We assume that H := £ {C*(n)}" and V £ 

Throughout this paper, we set x = {xi, Xn) G M", dj = for j = 1, 2, n 

and dt = ■§^, = (di, A = and JI ■ J denotes the scalar product of 

H,J £ R’". 

The main problems in this paper are 

Inverse source problem 

Let iJ, V, R,T C dil, T > 0 be given suitably. Determine f{x), x £ il. from y|rx(o.T)- 
Moreover we consider 

dtu{x,t) + H{x) ■ Vu{x,t) + V{x)u{x,t) = 0, x£il,0<t<T (1.3) 


and 


u{x,0) = a{x), X £ V,. (1-4) 

Inverse coefficient problem 

Let a and H be suitably given. Determine V{x) and/or Hix) by data w|rx(o,T)- 

Equations (1.1) and (1.3) are transport equations and are models in physical 
phenomena such as Liouville equation and the mass conservation law. Moreover the 
transport equation is related to the integral geometry (e.g., Amirov [1]). As for other 
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physical backgrounds such as neutron transport and medical tomography, see e.g., 
Case and Zweifel [ 8 ], Ren, Bal and Hielscher [18]. 

Our inverse problem is formulated with a single measurement, and Gaitan and 
Ouzzane [9], Klibanov and Pamyatnykh [15], Machida and Yamamoto [17] discuss the 
uniqueness and the stability for inverse problems for initial/boundary value problems 
for transport equations by Carleman estimates. The papers [15] and [17] discuss 
transport equations with integral terms where solutions y and u depend also on the 
velocity as well as the location x and the time t. 

The main methodology in [9], [15], [17] is based on Bukhgeim and Klibanov 
[7]. After that, there have been many works. Limited to hyperbolic and parabolic 
equations, we can refer for example to Baudouin, de Buhan and Ervedoza [3], Beilina 
and Klibanov [4], Bellassoued and Yamamoto [ 6 ], Imanuvilov and Yamamoto [11], 
[12], Klibanov [14], Yamamoto [21] and the references therein. Here we do not intend 
to give any complete lists of the references. In [9] and [15], the key Carleman estimate 
is the same as the Carleman estimate for a second-order hyperbolic equation and in 
order to apply the Carleman estimate one has to extend the solutions y and n to (1.1) 
and (1.3) to the time interval (—T, 0). Such an extension argument makes the proofs 
longer, and requires an extra condition to unknown coefficients and initial value as in 

[15] . In Sections 2 and 4, we prove Carleman estimates (Lemmata 1 and 3), which can 
directly estimate initial values. Thanks to our Carleman estimates, we can simplify 
the proofs of the stability and relax the constraints of the principal coefficients iL’s. 

As for inverse problems for transport equations with many measurements, see 
surveys Bal [2], Stefanov [20] and the references therein. Klibanov and Yamamoto 

[16] established the exact controllability for the transport equations by a Carleman 
estimate. 

Unlike [2] and [20], we discuss the inverse problems for a single initial/boundary 
value problem where we need not change initial values or boundary values. More 
precisely, in the formulation for the inverse problems in [ 2 ] and [ 20 ], we have to change 
boundary inputs on some subboundary and repeat measurements of the corresponding 
boundary outputs on other subboundary. One can apply the method of characteristics 
to the same kind of inverse problem for the first-order equation and see Belinskij [5], 
Chapter 5 of Romanov [19] for example. 

We set 


Q = O X (0,T) 


and 

f cIH-i- = {a; G dfl] {v{x) ■ H{x)) > 0}, 

\ = {a; € dfl] {i'{x) ■ h\x)) < 0}. 

Throughout this paper, we assume that ■0 G C^(fl) and If = (hi,...,hn) G 
{Ci(n)}” satisfy 

y := miii(If(x) ■ Vi/^(x)) > 0. (1.5) 

We here note by (1.5) that \H{x)\ 7 ^ 0 for a; G H. 

We give four cases where (1.5) holds. 

Case 1. We assume 


|Vd| >0 on H, H{x) = \7d{x) 


X G fl 
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with some d G Then (1.5) holds if we choose = d(x), x G il. 

Case 2. Let us assume that {{hi{x), ...,hn{x)); a; £ 11} C M" is separated from 
(0,...,0) by a hyperplane aiXi + ■■■ + anXn = 0 with some G M and |ai| + 

■ ■ ■ + \o.n\ 7 ^ 0. Then '(/’(a;) = aia;i H-+ a„a;„ or '0(a;) = —aia;i — • • • — a„a;„ satisfies 

(1.5). In particular, (1.5) holds if H{x) is a constant vector because {H{x)-, H} is 
composed of one point. In fact, the separation condition means that | [H (a;)-VV'(a;)) | = 

|aihi(x) H-+ anhn{x)\ > 0 for all x £ H or < 0 for all x £ 12. 

Case 3. Let 0 £ 12. We assume that there exists a constant (5o > 0 such that 


\H{x)\ > So, X £ 12. 

Then ipi^) = satisfies (1.5) if max^^.^^ |x| is sufficiently small. 

Proof. By the Cauchy-Schwarz inequality, we have 

n n n 

{H{x) ■ VV^(x)) = ^/if(x)^ + ^ he{x)'^xjdehj{x) 

j=i 

_i 

2 


( n / n \ / ^ 

fci \j=i ) 

>(5g -max|x|||iJ||{Loo(n)}n||ViL||{Loo(n)}nxn. 
xeo 

Therefore (1.5) holds true if 



h ( 


>min |iL(x)p — 

aiGr2 




i=\ 


X 

l=\ 


'xjdihj{x) 


a=i 


0 < m^ |x| < 




Case 4. We assume that there exists io £ {1, 2,..., nj such that hig (x) > 0 for all x £ 
12. Then we choose sufficienly small 6 £ K such that 12 C {(xi, X 2 ,..., x„); Xig > b}. 
Setting ^(x) = {xig — b)^, we verify that (1.5) holds. In fact, 


{H{x) ■ VV'(a;)) = 2hig{x){xig — b) > 0 for x £ 12 . 


Now we state the first main result concerning the stability for the inverse source 
problem. 

Theorem 1 

Let y G H^{Q) satisfy (1.1) and (1.2), and let (1.5) be satisfied with some constant 
fi > 0. We assume that 


i?(x, 0) 7 ^ 0, X £ 12 


( 1 . 6 ) 


and 


dty,dtRGH\Q), dtRG L^{0,T;L°^{n)). 
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Let 


T> 




(1.7) 


(i) We assume 


WdtvWmQ) < M 

with fixed constant M > 0. Then there exist constants 9 G (0,1) and C > 0, which 
are dependent on 0,r, iJ, M, \\dtR\\L^o,T-L^in)), such that 

{ e 

if f {H ■ u)\dtyfdS,dt] +1 f f iH-u)\dtyfdS,dt 
\Jo y \Jo JdQ+ 

for all / G L'^{n). 

(ii) Without the assumption in (i) concerning \\dty\\L^(^Q), there exists a constant 
C > 0, which is dependent on \\dtR\\L^o,T-L^in)), such that 

II/IIl^(0) < ( r [ \{H-i.)\\dty\^dS,dt 

\Jo Jon 




for all / G L^{n). 

(iii) In addition to (1.1) and (1.2), we assume 

y = 0 ondn_x{0,T). (1.8) 

Then there exists a constant C > 0, which is dependent on n,T,H, ||I^||L=o(n), 
11 ^^ 11 ^ 2 ( 0 ,T;L“(n)): such that 

C-i { r [ {H- v)\dty\^dS,dt\ < ||/|U2(n) <c{ T [ {H ■ ,.)\dty\^dS,dt 

\Jo JaQ+ J \Jo Jan+ 

(1.9) 

for all / G L^{n). 



The conclusion of (i) is a stability estimate of Holder type and holds under a 
priori boundedness ||9t?/llL2(Q) < M, which is called conditional stability. On the 
other hand, the conclusions of (ii) and (iii) are Lipschitz stability and in particular, 
with (1.8) we can have both-sided estimate (1.9) for our inverse problem. 

We apply Theorem 1 to the inverse coefficient problem of determining V{x). 
Theorem 2 

For j = 1, 2, let Uj,dtUj G H^{Q) and let 


dtUj + H{x) ■ Vuj + Vj{x)uj{x,t) = 0 in Q, (1-10) 


Uj{x, 0 ) = a(x) 


X G ft 


( 1 . 11 ) 
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and 



with suitably given a and h. We assume that there exists ib € satisfying (1.5) 

for H, (1.7) holds and 


Uj,dtuj e H\Q)nL^{ 0 ,T-,L^{n)), j = 1 , 2 . 


Moreover we assume 


II^IIl“( 0): l|5tWjl|L2(o,T;L“(n)) < M, j — 1,2 


(1.13) 


and 


|a| >0 on Q 


(1.14) 


where M > 0 is arbitrarily fixed constant. 

Then there exists a constant C > 0 depending on T, H, xp, a, M such that 


C ^||9t(Mi - u 2 )||L 2 (an+x(o,T)) < Wi - 1^2||L2(n) < C\\dt{ui - u 2 )||L 2 (an+x(o.T))- 

(1.15) 


Similarly to Theorem 2, we can discuss the determination of H, but we have to 
determine n functions as the components of H, and so repeats of measurements of 
boundary data after changing initial values suitably are necessary but arguments can 
be repeated similarly to the proof of Theorem 2. Here we omit detailed discussions 
for the determination of all the components of H, but we consider the determination 
of the potential in the case of potential flows. That is, we consider 


dtp{x, t) + div {p\7d{x)) = 0 in Q, 
p{x, 0) = a{x), X € ft. 


(1.16) 


The first-order partial differential equation in (1.16) describes the mass conservation 
under stationary potential flow V(i(a;). Then we discuss the inverse problem of deter¬ 
mining the stationary potential d. 

For the statement of the main result, we define an admissible set of unknown 
potentials d’s. Let constants M > 0, Sq > 0 and functions gi £ C‘^{dn), 52 £ (7^(911) 
be arbitrarily chosen. We define the admissible set of d’s by 


V = V{So,M,gi,g 2 ) := {d£ (^^(Q); |Vd| > 60 > Oon ft, 


IMIIc^(f2) — ^\dQ — 5ii 9[/d|gn — 52}- 


(1.17) 


For the function 32 given in (1.17), we set 


F+ = {a; £ 9H; g 2 {x) > 0}, 
F_ = {a; £ 9H; g 2 {x) < 0}. 


(1.18) 


We state our final main result. 

Theorem 3 

For j = 1, 2, let pj,dtPj £ H^{Q) and let 


dtPj {x, t) + div {pjVdj (x)) =0 in Q 


(1.19) 
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and 


Pj = h on r_ X (0,T), pj{x,0) = a{x), x G 

with suitable a and h. We assume (1-14), 

suPdgi,(max^gj 2 d(a:) - min^gj=j d(a;)) 


T> 


SI 


and 


dtPj G H\Q), 

\\Pj\\L‘^{o,T-,min)) < M, j = l,2. 


( 1 . 20 ) 


( 1 . 21 ) 


( 1 . 22 ) 


Then there exists a constant C > 0 depending on n,T,a,6o, M, gi,g 2 such that 
C~^\\dt{pi- P2 )\\l^(Tj,x(o,t)) < \\di - d2\\H‘^(n) < <^11194(^1 - p2)llL2(r_^x(o,T)) (1-23) 
for all di, ^2 G V. 


The proofs of Theorems 1-3 are based on an argument by the Carleman estimates, 
which was originated by Bukhgeim and Klibanov [7]. Here we used a modified ar¬ 
gument by Imanuvilov and Yamamoto [11], [12] which discussed for inverse problems 
for second-order hyperbolic equations. 

The paper is composed of four sections and an appendix. In Section 2, we prove 
a relevant Carleman estimate for the proofs of Theorems 1 and 2, and an energy 
estimate, and in Section 3, the proofs of Theorems 1 and 2 are completed. In Section 
4, we prove another Carleman estimate for the proof of Theorem 3 whose weight 
function gives a Carleman estimate also for the Laplacian and complete the proof of 
Theorem 3. In Appendix, we prove the Carleman estimate for the Laplacian. 

2. Key Carleman estimate and energy estimate. We recall that 

Q = H X (0,T) 


and we set 

Pu = dtu + H(x) -Vu + V{x)u, Pqu = dtu + H{x) ■ Vu, {x, t) £ Q, 


Mo = /3|[divi7||i«=(a)|[div {H{H ■ VV'))||L°°(n) 

and 

ip{x,t) =-f3t + i^ix), {x,t)GQ (2.1) 

with tp G C^(fi) and /3 > 0, and 

B(x) :=dt(p+(II-V(p) = -/^ + (II(x)-Vpj), xGfl. (2.2) 

First we prove 

Lemma 1 
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(i) We have 

Jn Jq 

<2 / + (sMq + 2||T7|||oo(n)) [ \u{x,t)\'^e^^‘^dxdt 

Jq Jq 

+s f [ B{x){v ■ H)\u\^e'^’^‘^dSxdt 

Jo Jan 

for all s > 0 and u € H^{Q) satisfying u{-,T) = 0 in fl. 

(ii) We assume (1.5) and 


0 < /3 < /r := min(iJ(x) • V'tp{x)). (2-3) 


Then 


s f \u{x,0)\'^e'^^‘^^^'°^dx + ^ f \u{x,t)\'^e'^^‘^dxdt 

Jn 2 Jq 


< 2 


\Pu\‘^e‘^^‘^dxdt 


s[ [ B{H-iy)\u\^e^^‘^dSxdt 

Jo Jdn+ 


for all s > So and u € ff^(Q) satisfying u{-,T) = 0 in il, where 


So = max 


4Mo A/8||f^|U°°(Q) \ 
[H-Py fi-P j' 


(2.4) 


Inequality (2.4) is an estimate of Carleman’s type, which holds uniformly for suffi¬ 
ciently large s > 0. We emphasize that Carleman estimate (2.4) can estimate also the 
initial value u(x, 0), and the weight function is linear in t. In [17] such a linear weight 
function is used for proving a stability estimate for an inverse problem for a transport 
equation. For the transport equation, the works [9] and [15] used the same weight 
function as the second-order hyperbolic equation, that is, ip{x,t) = gj'{d.{x)- 0 t )_ 
choice (2.1) of the weight function enables us not to need to take any extensions of u 
to (0, —T) in discussing the inverse problem. In [15], the extension argument requires 
an extra assumption on the initial value and unknown coefficients in addition to (1.14). 


Proof of Lemma 1 (i). 

First we assume that V = 0. We set w{x,t) = e'^‘^^^’*'>u{x,t) and {Lw){x,t) = 
e"‘^("^’*)Po(e“®‘^w). Then 


Lw = {dtw + {H{x) ■ Vw)} — sB{x)w. 
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Hence by u{-, T) = 0, we have 


f \Pou\^e^'^‘^dxdt = [ \Lw\^dxdt 

Jq jq 

= Idtw + {H ■'Vw)\^dxdt + / \sB\^w^dxdt — 2s / Bw{dtw + {H ■'Vw))dxdt 
J Q J Q J Q 

>—2s / B{dt'w + H ■'Vw)wdxdt + / B^w'^dxdt 

Jq Jq 

=—s I [Bdt{w'^) + BH ■y{w^))dxdt + j B'^w^dxdt 

Jq Jq 

=s ( {dtB + {diy BH))w^dxdt — s ( f B{v ■ H)w^dSxdt 
Jq Jo Jon 

+s^ / B^w^dxdt + s / B{x)\w(x,0)\^dx 
Jq Jn 

>—Mos [ w^dxdt — s [ f B{v ■ H)w‘^dSxdt 
Jq Jo Jan 

+s^ / B^w'^dxdt + s / B{x)\w{x,0)\^dx. 

Jq Jn 

Substituting w = e^‘^u, we have 

s f H(a:)|u(x, 0 )| 2 e 2 *^("’°)dx + s2 f B^{x)\u\'^dxdt 

Jn Jq 

< [ |PoMpe^®‘^c?xc?t + Mqs f \u\'^e^^‘^dxdt + s ( f B{x)\{H ■ v)\\u\^e^‘“^dSxdt. 
Jq Jq Jo Jan 

(2.5) 

Next let V € 0. Then 

|Pom|2 = \p^u + Vu- Fup < 2|Pou + Vu\^ + 2\Vu\^ 

<2|Pup + 2||y||i„(^)H^ 

Therefore (2.5) completes the proof of Lemma 1 (i). 

Proof of Lemma 1 (ii). 

By (2.2) and (2.3), we have B(x) > p, — /3 > 0 on O. Therefore we absorb the second 
term on the right-hand side of the conclusion of Lemma 1 into the left-hand side. 
More precisely, we have 

s^ [ B^{x)\ufe^‘‘'^dxdt>{n-I3fs^ [ \u\^e^^‘^dxdt, 

Jq Jq 


and so we easily see that if s > sq, then 
{d- P? 2 , ( (M-/3)\2 


>- 




s — Mqs ) + 


(M-/3)^2 


— 2||P||ioo(Q)^ 


2 
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Thus, by noting 


/an 


B{x){H ■ iy)\u\^dSxdt < 




B{x){H ■ v)\u\^dSxdt, 


the proof of the part (ii) is completed. 


Next we show the classical energy estimate. 

Lemma 2. 

Let |ll^||L~(n) < M and ||.ffll{c'i(o)}™ < M with arbitrarily fixed constant M > 0. 
Let w € H^{Q) satisfy 


f dtw + H{x) ■ Ww + Vw = F{x,t) in Q, 
[ ■u;(a;, 0) = a(a;), x G il. 


Then there exists a constant C > 0, depending on fl, T, M, such that 





{H ■ iy)\w\'^dSxdt 


( 2 . 6 ) 


<C\\\a 


lli^(n) + ll^lli=(Q)+^^^^ \{H-v)\\w\^dSxd^, 0<t<T. (2.7) 


Proof. 

Multiplying dtw + H- 'Vw + Vw = F hy 2w and integrating over fl, we have 
dt f \w{x,t)\'^dx + '^ f hidi{\w\'^)dx + 2 ( V\w\^dx = 2 f Fwdx. 

Jq Jq Jq 

Setting E{t) = \w{x,t)\'^dx and integrating by parts, we obtain 

E'{t) = — f {H ■ i>)\'w\‘^dSx + f {div H)\'w\^dx 
JdQ Jn 

—2 / Vw'^dx + 2 / Fwdx. 

Jq Jq 

Therefore, noting that 2 \Fw\dx < \F\’^dx + \w\'^dx and integrating over (0, t), 

we have 

E{t)-E{0) = - [If + f ] {H-iy)\wfdSxdt+ [ f {div H)\w\‘^ dxdt 
Jo \Jdn+ JdQ-J Jo Jn 

—2 [ [ Vw^dxdtF [ [ \Ffdxdt+ [ E{^)d^. 

Jo Jn Jo Jn Jo 


Therefore 


-E'(^) + j < ||a||^2(Q) + ||F||i2(Q) 
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+ (2||y|Uo.(0) + ||diviJ|U^(0) + l) / E{C)d^ 


{H ■ i')\w\^dSxdt 


0 JdQ.- 


( 2 . 8 ) 


for 0 < t < T. Since f* fgQ^{H ■ i')\w\‘^dSxdt > 0, 0 < t < T, omitting the second 
term on the left-hand side of (2.8), we apply the Gronwall inequality, and we obtain 


E{t) < ^||a|||2(Q) -I- ||^|li2(Q)) 


+ [ [ \{H ■ i^)\\w\^dSxdt\ 0<t<T. (2.9) 

Jo JdQ- ) 

Substituting (2.9) into the third term on the right-hand side of (2.8), we complete the 
proof of the lemma. 

3. Proofs of Theorems 1 and 2. 

Proof of Theorem 1. 

Henceforth C > 0 denotes generic constants which are independent of s > 0. We 
note that Py = dty -f El ■ Vy + Vy and y, = min^gj2(i7(a;) • V^(x)) > 0, and we set 
R = max^gQ ^^(a;) and r = min^g^'(/;(x). By (1.7), we can choose /3 > 0 such that 

T>^^, 0</3<m. (3.1) 

In fact, it suffices to choose /3 > 0 such that /3 G (0, y) is sufficiently close to y. 

With this /3 > 0, we set 

(p(x,t) =—/3t + (x,t)GQ. (3.2) 


Then (3.1) implies 

ip(x, T) < R — [3T < r < (p{x', 0), x, x' G fl. 

Therefore by y? € C^{Q), there exist (5i > 0 and tq, ri such that R —(3T < ro < ri < r. 


f y}[x,t) > ri, X G n, 0 < t < i5i, 

\ V?(x, t) < tq, X G El, T — 25i <t<T. 


For applying Lemma 1, we need a cut-off function y G C^(R) such that 0 < y < 1 
and 



0 < t <T-26i, 
T - 6i <t <T. 


(3.4) 


We set 


2 = {dty)x- 


Then z{x, T) = 0, x G El and 


Pz = xfidtR) + {dtx)dty, (x, t) gQ 
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and 


z{x,0) = f{x)R{x,0), X € ft. 

Applying Lemma 1 (ii) to z, we obtain 

s[ \z{x,0)\^e^^‘^^^’°Ux<C [ \xf{dtR)\^e^‘’^dxdt+C [ \{dtx)dty\'^e'^^‘^dxdt+Ce^^D^ 

Jn Jq Jq 

(3.5) 

for all large s > 0. Here 

D^= [ [ {H ■iy)\dty\^dS^dt. 

Jo J 

Since dtX = 0for0<t<T — 2Si or T — Si < t < T, hy (3.3) and the a priori bound 
WdtvWmQ) < M, we have 


JQ 

and 


\{dtx)dty\^e^'"'^dxdt = 


f-TSi 


f T —2^1 J O 


\{dtx)dty\^e^‘^^dxdt < 


r-T-5i 


' T — 26i J 


\{dtx)dty?e^^'f^dxdt < 


\dty\'^dxdt 
(3.6) 


(3.7) 


for all large s > 0. Moreover R{x, 0) yf 0 for x £ fl and z{x, 0) = f{x)R{x, 0), x £ fl, 
we have 

[ \z{x,0)\^e'^^‘^^^’°Ux>C [ \f{x)\^e'^^'^^^’°Ux. 

t/ O Q 

Therefore (3.5) yields 

s[ |/(a;)|2e2*‘^(^’°)dx < C / |/(x)pe2'*^(“’‘)dxdt + 

Jo Jq 


Since yy{x,t) < </3(a:, 0) for {x,t) £ Q, we have 

/o Jn 


s[ |/(a;)|2e2*‘^(^>°)dx < C / [ \f{x)\‘^e^^‘^^^’°^dxdt + CM‘^e‘^^^^+Ce^W^ 

</ Q J 0 J ^ 

<CT [ \f{x)\^e^‘“^^^’°'^dx + CM^e^^^°+Ce^W^, 

Jq 

that is, 

{s-CT) [ |/(a;)pe2®‘^(“’°)da; < 

Jq 

for all large s > 0. Using ^p{x, 0) > ri by (3.3) and choosing s > 0 large, we obtain 

[ \f{x)\^dx<CM^e^^^°+Ce^^D'^, 

Jq 
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that is, 


ll/lli=(n) < (3.8) 

for all large s > where s* > 0 is a sufficiently large constant. Here we set 
r* := ri — ro > 0. We separately consider the two cases: D > M and D < M. 

Case 1 D > M: 

Estimate (3.8) implies 

ll/lli=(0)<(^^e-2-+Ce<^«)i5^ (3.9) 

Case 2 D < M: 

Replacing C by , we see that (3.8) holds for all s > 0. We make the right-hand 

side of (3.8) small in s. We choose that is, 

2 , M 


Therefore (3.8) is reduced to 


||/|U2(n) < 2CM^-^dO, 


where 9 = ^ Thus the proof of Theorem 1 (i) is completed. 

Next we prove Theorem 1 (ii). We have 


f dt{dty) + H{x)-Vdty+ Vdty = f{x)dtR, {x,t) G Q, 
1 {dty){x, 0) = f{x)R{x, 0), x G H. 


Applying Lemma 2 to dty, we obtain 


\dty{x,t)\^dx 


/ 0Q4 


{H ■ v)\dty\^dSxdt 


<C\\f\\l.(^^)+C [ [ \{H ■ i.)\\dty\^dS.dt (3.10) 

Jo JdQ- 

for 0 < t < T. Therefore, omitting the second term on the left-hand side of (3.10) 
and applying it to (3.6), we obtain 


/ |(5a)9t2/pe2®‘^dxdt<Ce2®’'«||/||i.(^)+Ce2®’'o / [ \{H ■ i.)\\dty\^dS, 

Jq Jo JdQ.- 

and similarly to (3.8), from (3.5) we can obtain 

^ [ \{H ■ i^)\\dty\'^dSxdt. 


dt 


\Un) < 


+ C'e 


Cs 


lo JdQ 


Choosing s > 0 large, we can absorb the first term on the right-hand side into the 
left-hand side, and complete the proof of (ii). 
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Finally we prove (iii). By (1.8), the conclusion of (ii) immediately yields the 
second inequality in (1.9). Next, by (1.8) and (3.10), we have 

Jo Jan+ 

which proves the first inequality in (1.9). Thus the proof of Theorem 1 (iii) is com¬ 
pleted. 

Proof of Theorem 2. 

Theorem 2 can be derived directly by Theorem 1. In fact, setting y = ui — U 2 -, 
f = Vi — V 2 and R = —M 2 , by (1.10) - (1.14) we have 

( dty + H{x)-Vy + Viy = f{x)R{x,t), {x,t)eQ, 

I y{x,0) = 0, X e fl, 

[ y|aa_x(o.T) = 0, 

and dty G \\dtR\\L'^(o,T-,L°°in)) < R{x,0) = -a{x) 0 for x G O. Thus 

Theorem 1 (iii) yields the conclusion of Theorem 2, and the proof of Theorem 2 is 
completed. 

4. Proof of Theorem 3. 

We set 


PdU = dtu + Vd ■ Vu + uAd. 


For the proof, we further need a Carleman estimate for A with the same weight 
function for Pd- Unfortunately the weight function defined by (2.1) does not work 
as weight for a Carleman estimate for A. Thus we have to introduce a second large 
parameter in the weight function. That is, as the weight function, we set 

(^d(x,f) = e^(-^‘+‘'("», (x,t)Gg, (4.1) 


where A > 0 is chosen later. The weight function in the form (4.1) has been known as 
more flexible weight function producing a Carleman estimate (e.g., Hormander [10], 
Isakov [13]). 

In the existing works on inverse problems by Carleman estimates, given weight 
functions have been used, and to the best knowledge of the authors. Theorem 3 is 
the hrst case where a Carleman estimate with unknown coefficient as the weight, is 
seriously involved. The use of such a Carleman estimate is necessary in order that the 
admissible set V of unknown coefficients d’s is generously formulated such as (1.17). 
Otherwise the admissible set of unknown coefficients should be more restrictive. For 
example, for inverse problems of determining principal parts for second-order hyper¬ 
bolic equations, there are no works by Carleman estimate with weight function given 
by unknown coefficient, so that choices of admissible sets of unknown coefficients of 
the principal parts are very limited (e.g., Bellassoued and Yamamoto [6], Section 6 
of Chapter 5). For it, we need a Carleman estimate where all the constants can be 
taken uniformly for arbitrary d G P. 

We first show a Carleman estimate for Pd with weight (4.1), which should hold 
uniformly for all d G 17. We fix /3, T > 0 such that 


0<P<6l 


rj. ^ sup<igp(max^gjj d(x) - min^gj^ d(x)) 


(4.2) 
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where the constant characterizes V (see (1.17)). Henceforth we set 
Jd{x, t) := dt<^d + (Vd • V(/9d) = Xipd{-I3 + |Vd|^). 


Lemma 3 

For each A > 0 there exists a constant sq > 0, which is dependent on X,So,M, such 
that we can choose a constant C > 0 satisfying 


sAi^d(x, 0)|it(a::, X‘^ip^\u{x,t)\'^e^^‘^‘^dxdt 

1 JQ 

<C [ \Pdu\^e^‘^^dxdt+ [ [ sJd{d^d)\u\^e^^‘^‘‘dS^dt (4.3) 

Jo Jo Jon 


for all s > So, all d € 27 and all u G H^{Q) satisfying u(x, T) = 0, x € H. 


Here we note that the choices of C and sq in (4.3) are uniformly for d £ 27. 

Proof. 

The proof is similar to Lemma 1 by noting the independency of the constant of d £ 27. 
Let d £ 27 be chosen arbitrarily. First we consider the case of PdU = dtu + Vd • Vu. 
We set w = and 


Lw = e®‘^'^Pd(e = dtw + Vd • Vw — sJd{x, t)w. 

We note dtJd{x,t) = -X^l3ipd{-P + |Vdp), dkJd{x,t) = X^ipd{dkd){-(3 + |Vdp) + 
A(^d(9fe|Vdp) and 

Jd(.x,t)>X^pd{5l-(3), {x,t) gQ, dGV (4.4) 


by (1.17). Then 

f |Pd'upe^'“^'^dxdt = f \Lwfdxdt 

Jq Jq 

= / \{dtw+ X/d-X/w) — sJdw\^dxdt 

Jq 

= / \dtw + X/d ■ X/wf dxdt + / J2\wfdxdt — 2s / {Jdw){dtw + X/d ■ X/w)dxdt 

J Q J Q J Q 

>s^ / dJiwpdccdt —s / Jd{2w{dtw))dxdt — s / JdX/d ■ {2wX/w)dxdt 

Jq Jq Jq 


='-Ii +12 + Is- 


Henceforth C > 0 denotes generic constants which are independent of d £ 27 and 
s > 0, but dependent on Sq,M,X. Next integration by parts yields 

l 2 = -s / Jddt{\w\^)dxdt =-s / [Jd|typ]o dx + s / (9* Jd)|ri;pda;dt 
J Q J Q, J Q 

>s j X(pdix,0){—l3 + \X7d\'^)\w{x,0)\'^dx — C j sX'^(pd\w\'^dxdt 
Jq. Jq 

>s((5q —/3) f sX(pd{x,0)\w{x,0)\'^dx — C ( sX^ipd\w\^dxdt 
Jq Jq 
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and 


„ n 

h = -s / y^{Jddkd)dki\w\'^)dxdt 
JQ k=i 


= — s 


/o JdQ 


'^Jd{dkd)iyk\wfdS^dt + s / dk{Jddkd)\'wfdxdt 


k=l 


Qk= 


k=l 


> - 


lo Jdn 


sJd{d^d)\'w\'^dSxdt — C / sX^ipd\'w\'^dxdt. 

Jo 


Hence 


sA(/5d(a;,0)|w(a;,0)fdx + / dxdt 


< C / \Lw\^dxdt + C f sX'^(fid\w\'^dxdt + ( f sJd{d,yd)\w\^dSxdt. (4.5) 

Jo Jo Jo JdO, 


We choose 

so > 

Then, since (pd{x,t) > for {x,t) € Q and d € V, if s > sq, then 

C 


s^X^^l - CsX^^d = s"AV^ ( 1 - 


2 \ 2,„2 


S(pd 


>s^X\l [ 1 - 


^gA(/3T+M)x 1 „ „ „ 

-^- > T^S^X^V’l 

So J 2 


Therefore (4.5) yields 


sXipd{x, 0)|r(;(x, 0)\'^dx 


s^ X^ ip^\wf dxdt 


f \Lw\'^dxdt+ f [ sJd{dvd)\w\^dSxdt (4-6) 

Jq Jo Jon 

for all s > So- Henceforth So > 0 denotes generic constants which are dependent on 
X,6o,M. Since w = we rewrite (4.6) in view of u, so that 

[ sX(pd{x,0)\u{x,0)\^e^^‘^‘‘^^’°'^dx + [ s^A^(/9^|Mpe^®‘^‘^(ixdt 


<C [ \Pdufe^^‘^-^dxdt+ [ [ sJd{d^d)\u\^e^^‘f’’‘dS, 

Jq Jo Jdo. 


dt 


for all s > Sq. Here we note that C > 0 and so > 0 may vary line by line, but 
dependent only on X,6q,M. Next 

\dtu + Vd • Vu + (Ad)Mp < 2\dtu + Vd • Vwp + 2|(Ad)up 
<2\dtu + Vd • Vup + 2M2 |m|2 
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in Q for each d G T). Similarly to the proof of Lemma 1 (i), we can finish the proof of 
(4.3). 


Next we show a Carleman estimate for A with the weight function ipj_{x,0). 

Lemma 4 

There exists a constant Aq = Ao(i5o, M) > 0 such that for any A > Aq, there exists a 
constant si = si(A, (5o,M) > 0 satisfying: we can choose C > 0 such that 

[ (sA^(pd(a:, 0)|V/(a:)p + s^A^<pd(a;, 

Jn 

<C [ 

Jq 

for all s > si, d G V, f G Hq{Q). 


Lemma 4 is a Carleman estimate for A and the Carleman estimate for A is well 
known (e.g., Hormander [10], Isakov [13]). However we need the uniformity of the 
constants so,C with respect to d € 2? in the Carleman estimate. Thus in Appendix 
we give a direct proof by integration by parts, which differs from [10], [13]. 

Now we proceed to 

Proof of Theorem 3. 

First Step. We set 

mg = sup(m^d(a;) — mind(a;)). 

dev xeQ xev 

By assumption (1.21): T > ^ of the theorem, for arbitrary d G V, we can choose 
/? > 0 and fix such that 

0<f3<6l > ^- ( 4 - 7 ) 

We choose Ag > 0 given in Lemma 4 and fix. Next we set s* = max{sg(Ag, dg, M), si(Ag, dg, M)}, 
where Sj(Ag, dg, M), j = 0,1 are proved to exist in Lemmata 3 and 4. Since d^d = g 2 
on dVl for any d G V, we note that the subboundaries T.). and r_ defined by (1.18) 
correspond to and dfl- respectively. Since Ag > 0 is fixed, henceforth by C 
denoting generic constants which are independent of s > s* and d G T) but dependent 
on Ag, a, /3, mg, O, T, noting (4.7) and (4.4), by Lemmata 3 and 4 we obtain 

f s\y{x,0)\'^e^^‘^'‘^^’°'^dx + f s^\y\^e^^‘^'‘dxdt 
Jn Jq 


<C \Pdyfe^^'^'^dxdt + Ce 


Cs 


'Q 


\y\^dS^dt 


0 


(4.8) 
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for all s > s^:, d € V, y G H^{Q) satisfying y{-,T) = 0 in fl and / € Thus we 

obtained two Carleman estimates which hold uniformly for arbitrary d G T>. 

Second Step. We set y = pi — p 2 , f = di — d 2 and R{x,t) = —p 2 {x,t) in Q. Then, 
by (1-19) and (1.20), we obtain 

r dty + Vdi • Vy + yAdi = V/ • Vi? + RAf in Q, 

< yix,0)=0, xG^, 

[ y = 0 on r_ X (0, T). 

We differentiate this partial differential equation with respect to y in t, noting that 
dty G H^(Q), dtR € L^(0,T;lT^’°“(n)) and R{x,0) = —a{x). Then we have 

f dt{dty) + Vdi • V(dty) + {dty)Adi 
I = Vf-VdtR+(dtR)Af inQ, 

) dty(x,0)=-Vf-Va-aAf, x G fi, ^ 

[ dty = 0 onr_x(0,T). 

Applying Lemma 2 to (4.10) and estimating dty, we readily verify the first inequality 
in (1.23). Thus the rest part of this section is devoted to the proof of the second 
inequality of (1.23). 

Since y{-,T) = 0 does not hold, for applying (4.8), we need a cut-off function. 
Noting 


and 

ipdx{x,T) = < g--^o/3TgAo max„,gjjdi(a:) 

for all a; € n. Henceforth we set 

ro = ^Aoe-^“'5-^«'^(/3r-mo)- 

By (4.7) we have rg > 0. We note that rg > 0 is independent of s > 0 and special 
choices of di, ^2 € 77, but dependent on Sg, M, Ag. 

The mean value theorem implies 

maXipdAx,T) - minipd,ix,0) = _ ^Ao 

=e^'^Ao{(m^di(x) — min(ii(a:)) — /3T} < Aoe^^(mo — /ST), 

xeQ x£Q. 

where is a number between Ag min^j.^^ di(a::) and —XgfdT + Ag max^^Qdi{x). 
Therefore, by d G 27, we have > g-^odT-XoM^ 

inax<Pdi{x,T)-min(pdAx,0) < -Age^^(/3T-mg) < {(dT-mg) = -2rg. 

(4.11) 

For any di G 27, we have 

I va.&xipdAx,T)—xaaxLpdAx,t)\ = \ exp(Ag(—/3T+m^di(a;)))—exp(Ao(—/3t-|-m^di(a;))| 

x^tl x^tl x^il x^Q, 


= I exp(Ag m^di(a;))||e — e < e'^“^Ag/3|r — t|. 


(4.12) 
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At the last inequality, by di € 2? and the mean value theorem, we have 

m‘Axdi{x) < M 

xGO, 


and 


|g-Ao^T _ g-Ao/3t| ^ |e-52(_Ao/3r + Ao/3t)| < Ao/3|T - t|, 

where ^2 > 0 is some constant. 

We fix a constant (5i > 0 sufficienly small such that 

2e^°^XoPSi < ro. 

Then, for T — 26i < t <T and any di G V, by (4.11) and (4.12), we have 

max ipdi (x, t) - min (x, 0) 
xen a;Gn 

= mEix (pdi (x, t) - max pdi (x, T) + max ipdi (x, T) - min pdi (x, 0) 

xEQ 

<1 maxv5di(x,t) - inaxipdiix,T)\ - 2ro < e^°^Xo(3\t - T\ - 2ro 
< 2 e^°^X 0 P 61 - 2 ro < cq - 2ro = -xq. 


Therefore 


maxv5di(x,t) < ^0 -'To, T-2Si<t<T, (4.13) 

where we set po := sup^^g^, min^go 

Let X G C^(]R) satisfy 0 < x < 1 and (3.4) with here chosen 5i. We set 

2 = xdty- 


Then, by (4.10), we have 

{ Pd^z = dtz + Vdi • Vz + zAdi 

= f ■'^9tR + x{dtR)Af + {dtx)dty in Q, 

z(x,0) =-V/-Va-aA/, z(x,r) = 0 x £ O, 
z = 0 onr_x(0,T). 

Applying (4.8) to (4.14), we obtain 


(4.14) 


[ s|aA/ + V/- 

JQ. 


<c [ \xyf-VdtR + x{dtR)Af\'^e^^’^'‘^dxdt + C [ \{dtx)dty\^e^^’^'‘^dxdt + 

JQ JQ 

(4.15) 


for all s > s*. Here and henceforth we set 

pT r 


= 


\dt{pi - p2)\ dSxdt. 


0 
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Apply Lemma 2 to (4.10), by (1-21) we have 

^ \dty{x, t)\'^dx < C(|| V/ • Va + aAf\?mn) + IIV/ • V(5*A) + {dtR)M\\h^Q)) + 


^ C'(||A/|||2(q) + ||V/|||^2(f2)}n + 0 < t < T. (4-16) 

Consequently, since dty ^ 0 only if T — 25i < t < T — 5i hy (3.4), it follows from 
(1.14) that (4.15) and (4.16) yield 


Jq Jq 

<C [ (|A/p + 

JQ 

\Jt-2Si Jq J 


for all s > s*. For the first integral on the right-hand side, we used (pd-^{x,t) < 
ipdi{x,0) for {x,t) £ Q, and so 


/ IxV/ • VidtR) + x{dtR)Af\^e^^'^‘‘^’^^’*'>dxdt 
Jq 

<C [ (|V/p + |A/ne2'*^‘^i(“’‘)rf2.rfi<C' [ (|V/p-L|A/ne2«‘^"i("^’°)da;. 

Jq Jq 

Therefore we obtain 

Jq 

<C [ Cs [ 

Jq Jq 

+Ce^^D^ 

for all s > s^,. The first term on the right-hand side can be absorbed into the left-hand 
side by choosing s > 0 large, and by (4.13) we have 


< Cs + C'e2«(/^o-’'o)(||A/||i2(n) + II V/||^^ 2 (n)}.) + 

(4.17) 

for all s > s*. 

Since / = di — ^2 £ ^o(^) by di,d 2 £ T), we can apply (4.9) to obtain 

[ (s^lv/p-f < C [ s\Af\^e^^‘^'‘^^^’°'>dx. 

Jq Jq 
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Substituting this into the left-hand side of (4.17), we obtain 

[ (s|A/|2 -f s^\Vf\^ + 

Jn 

for all s > s*, and absorbing the first term on the right-hand side into the left-hand 
side again by choosing s > 0 large if necessary, we have 

Se^®'^“(||A/||^2(Q) -I- ||V/|||^2(Q)}n -I- ||/||i2(Q)) 

<s [ (|A/|2-F|V/p-f 
Jn 

<^g2.(^o-.o)(|| 

for all s > s*. Noting C'-i||/||^ 2 (n) < ||A/||| 2 (n) + IIV/||^i 2 (n)}n + Il/Ili 2 (n) < 
C||/llff 2 (Q), we see that 

||/lk^(0) < 11/11^2(0) 

for all large s > s*. Since ro > 0, for large s > 0, we can absorb the first term on the 
right-hand side into the left -hand side, so that 

ll/ll«^(a) 

Thus the proof of Theorem 3 is completed. 


Appendix. Proof of Lemma 4 

In order to clarify the dependence in the Carleman estimate on a-nd M, for arbitrary 
d € T>, we prove the lemma directly and the proof is similar for example to Section 3 
in Yamamoto [21], where a Carleman estimate is proved for a parabolic equation. 

By the usual density argument, it suffices to prove for / € Cq)!!). We set 

ipd{x) = ifd{x, 0), g{x) = <7 = A/ 


and 


Log{x) = 


We have 

Loff = Ag — 2 sA^(iV(i • V 5 -t-s^A^^^lVdpg — sA^V'dlVdpg — sA' 0 d(Ad )5 in fl. ( 1 ) 
We set 

Ai := -sAVdiVdp - sXMM- (2) 

We note that 

|Ai(x, s, A)| < CsX'^tpd for a; G fl and all large A > 0 and s > 0. 
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Here and henceforth by C, Ci , etc., we denote generic constants which are independent 
of s, A, tpd but dependent on M, Jq > 0, and these constants may change line by line. 
Moreover we always assume that A > 1. Hence we use the following inequality: 
A™ < A”"' if 0 < m < m'. 

In view of Hi, we have 

Log = Ag-2sAV^dVd-Vg + s^AVd|Vd|^g + Hig = ge'*^‘' in H. (3) 


Taking into consideration the orders of (s, A,^c;), we split L^g as follows: 


I Lig = Ag + s'^X‘^ip'^\Vd\‘^g + Aig, 

\ L2g = — 2 sA' 0 (/V(i • Vg in H. 

Since = \\Lig + L2g\\l2(^a)^ obtain 

2 f {Lig){L2g)dx < f q^e^’^'^'^dx. 

JQ JQ 

Now we estimate 


/ {Lig){L2g)dx 
Jn 

= - f 2sAV^d(Vd • Vg)Agdx - [ 2s^X^iljl\\7d\^{\7d ■ \7g)gdx 
J Q J Q 


(4) 


(5) 


3 

— / 2sXipdAi{'X/d ■ 'X/g)gdx =: 'S^ Jfc. 

Jn 


By the integration by parts and g = d^g = 0 on dXl, we have 


( 6 ) 


Ji = / 2sA'0d ^ {dkd){djdkg)djgdx + / 2sA ^ dj{{dkd)ipd)idjg)dkgdx 


jA=i 


jA=i 


= / sXipd E {dkd)dk{\djg\‘^)dx + / 2sA ^ dj{{dkd)ipd){djg)dkgdx 


iA=i 




P n p n 

/ sX^dk{{dkd)il)d)\yg\’^dx + / 2sA ^ dj{{dkd)'iljd){djg)dkgdx. 

fe=i 7,fe=l 


Since 


dj{{dkd)'tpd) = {djdkd)'ijjd + {dkd)djtpd = X{djd){dkd)iljd + {djdkd)'ijjd, 
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we have 


Ji = — [ sX^ilJd\yd\'^\Vg\'^dx — f sX{Ad)'ipd\'^g\^dx 

Jn Jn 

P n p n 

+ / 2sX‘^'ipd {djd){dkd){djg)dkgdx + / 2sXipd {djdkd){djg)dkgdx 

Jn Jn 

= — f sX^'tpd\'^d\'^\Vg\'^dx — f sX{Ad)'tlJd\'^g\^dx 
*/ r2 r2 




2sX^'ipd 


J2 ^^3d)djg 


j,k=l 


dx+ 2sX'ipd E {djdkd){djg)dkgdx 




>— f sX'^ipd\'^d\'^\Vg\'^dx — j sX{Ad)'tpd\'^g\^dx 
Jn Jn 


2sX-ipd y] {djdkd){djg)dkgdx. 


j,fe=i 


Therefore 


Ji> — [ sX'^'tpd\'^d\'^\Vg\'^dx — C j sA'i/'dlVgpdx. 
Jn Jn 

Moreover we have 

P n 

J2 = - / s^X^'ipl\X/d\'^'^{dkd)dki\g\‘^)dx 

k=i 

P n p n 


(7) 


fc=l 


/c=l 


and 


= f 3s^AVd|Vd|^|5pdx-C / s^X^'tpl\g\'^dx 

J n J Cl 

P n 

J 3 = - / sXipdAi y^^{dkd)dki\gf)dx 

k=i 

P n p n 

= / '^sXdk{,'>l3dAi){dkd)\g\^dx + / sXil)dAi{dld)\g\^dx 


‘ /c=l 


>-c[ s^X'^ijl\g\^dx. 

Jn 

Therefore ( 6 ) - (9) yield 

\ [ > f {Lig){L2g)dx 

^ Jn Jn 

> — / sX'^'ipd\^d\‘^\X/g\'^dx + f 3s^X‘^ipd\^d\‘^\g\'^dx 
Jci J Cl 


( 8 ) 


(9) 
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-C j s\il)dWg\^dx-C [ +s’^X^'ipDlgl'^dx. (10) 

J Q J Q 

For the proof of the lemma, we have to estimate sX^'tpdl'^g]^ + s^X‘^'ip^\g\'^, but the 
first and the second terms on the right-hand side have different signs. Thus we need 
another estimate. That is, multiplying (3) with —sX^ipdl^dl^g and integrating over 
n, we obtain 

^ 4 :=- [ {Ag)sX'^'ilJd\yd\'^gdx+ [ 2s'^X^iljl\yd\'^{yd ■ Vg)gdx 

j _O <^0 

- [ s^X‘^ipd\'^d\‘^\g\‘^dx - [ sX‘^ipdAi\X/d\‘^\g\‘^dx 
Jn JQ, 


= — sX^-ipdl^dl^ gdx. 

Jn 


( 11 ) 


Hence, by integration by parts, we obtain 

Ii= [ sX^^d\yd\^\yg\^dx+ f sA^VdVdlVd) • (V5)ffdx. 

Jn Jn 

Here we have 

[ sX‘^X/{\X/d\‘^tlJd) ■ ig'^g)dx < C [ sX^'tpdlgW'^gldx 
Jn Jn 

=C [ (sAVd|5|)(A|V3|)dx < I [ {s^X'^^llgf + X^\ygf)dx. 

Jn ^ Jn 

Consequently 

h > f sA^^/’d|Vdp|V 5 pc?a; — C f s^A^^/’^|gpc?a; — C f X'^\X/g\'^dx. (12) 

t/ Q t/ Q fi 


Next 


72 = / s^X^^ljl\VdWx/d-V{\gf))dx 
Jn 

« n „ n 

= - / s"A3^afc(V-^)|Vd|2(5fcd)|g|2dx- / s2AVdI]5fc((5ferf)|Vd|")|<?|"dx 


fc = l 




>-cf {s^X'^il^l X- s^X^il}l)\g\^dx >-C f s'^X*ij'^\g\^dx. (13) 

Jn Jn 


By (2), we see 


/i > —C j s^X'^'ijj‘]\g\^dx. 

Jn 


(14) 


[ qe^^'^sX^U'^dl^gdx <]- [ \q\^e^-'>’^dx+\ [ X*^|;l\yd\^\g\^dx, 

Jn J Jn J Jn 


Since 
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it follows from (11) - (14) that 

f — f s^X'^'il}^\Vd\^\g\'^dx — C f X^'il}^\g\^dx — C f X^\yg\^dx 

«/0 «/f 2 


/n 

Therefore 


f sA^'i/'d|Vdp|V5pdx — f s^X'^'tlj^\Vd\'^\g\'^dx 
Jq Jq 

<i / + C f s^X^ip^\g\^dx + C f X^\X/g\^dx. (15) 

2 Jo Jo Jo 

Thus we consider (10) + 2x (15): 

f sA^^dlVdplV^pda: + f s^X'^ip^\X/d\^\g\'^dx 
J Q J Q 


-^ [ + C / (sAi/;(i + A^)|V 5 pd 2 ; + C [ {s^X^tj}d +s^^^'4’d)\9\^dx. (16) 

2 Jo Jo Jo 


Since IV'jCj^)! > e for all a; G fl and d G we have 


Cni,-^ r<p>^M 

CX^ = sX^ipd X -^ < sX'^ipd -; 


(17) 


CsXtIJd = sX^ipd^, Cs^X^ipl = s^X'^'ipl^ (18) 

and 

c»nv3 = sH'i,i 9 *£- < .nv-3—^ (19) 

s s 

We choose sufficiently large Aq = Xq{5q,M) > 0 such that |^ < |; in (18). For any 
given A > Aq, we choose si = si(A,(5o,M) such that — < A in (17) and (19). Then 

CA^ CsXi^d < ^sAVj, 

Cs^X^iIjI, Cs^X'^tpj < 

for all A > Ao and s > si, and noting |Vd| > do > 0 for all d G TX, we can absorb the 
second and the third terms on the right-hand side of (16) into the left-hand side to 
obtain 

1 [ sX‘^-ipd\^g\^dx+^ [ s^X‘^'ip^\g\^dx<C [ ( 20 ) 

2 Jo 2 Jq Jo 

We re-write by means of /. Replacing g = we have jgp = and 

|y/|2g2«V-. = |y^ _ sA(Vd)V^de*’^Vl' 

<2|V5p -f 2s2A2|Vdp,/i2|/|2e2'*bd < 2|Vgp + Cs^A^dI 

Substituting them into (20), we complete the proof of Lemma 4. 
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